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Abstract 

CO A linear singularly perturbed convection-diffusion problem with characteristic 

layers is considered in three dimensions. Sharp bounds for the associated Green's 
function and its derivatives are established in the Li norm. The dependence of 
these bounds on the small perturbation parameter is shown explicitly. The obtained 
estimates will be used in a forthcoming numerical analysis of the considered problem. 
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^ 1 Introduction 

^ In this paper we consider the following problem posed in the unit-cube domain Vl = (0, 1)'^: 

2 Lxu{x) = ~e Axu{x) — dx^{a{x)u{x)) + b{x)u{x) = f{x) ioi x e il, (1.1a) 

^ u{x) = ioTxedn. (1.1b) 

>- Here e is a small positive parameter, and we assume that the coefficients a and b are 

^ sufficiently smooth (a, b e C°°{fl)). We also assume, for some positive constant a, that 

^ a{x) > a > 0, b{x) - 9^,a(x) > for all x e n. (1.2) 



Under these assumptions, (1.1a) is a singularly perturbed elliptic equation, also referred to 



as a convection-dominated convection-diffusion equation. Its solutions typically exhibits 
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Figure 1: Anisotropy of the Green's function G associated with (1.1) for e 
X - 

two-dimensional graph for fixed ^3 = X3 



0.01 and 



|, ^, |). Left: isosurfaces at values of 1, 4, 8, 16, 32, 64, 128, and 256. Right: a 



sharp interior and boundary layers. This equation serves as a model for Navier-Stokes 
equations at large Reynolds numbers or (in the linearised case) of Oseen equations and 
provides an excellent paradigm for numerical techniques in the computational fluid dy- 
namics 



The Green's function for the convection-diffusion problem (1.1) exhibits a strong 



anisotropic structure, which is demonstrated by Figure [Tj This reflects the complex- 
ity of solutions of this problem; it should be noted that problems of this type require 
an intricate asymptotic analysis [121 Section IV. 1], [13]; see also [201 Chapter IV], [TH 
Chapter III.l] and [HI [T3]. We also refer the reader to Dorfier [Ij, who, for a similar 
problem, gives extensive a priori solution estimates. 



Our interest in considering the Green's function of problem (1.1) and estimating its 



derivatives is motivated by the numerical analysis of this computationally challenging 
problem. More specifically, we shall use the obtained estimates in the forthcoming paper 
[0] to derive robust a posteriori error bounds for computed solutions of this problem using 
finite-difference methods. (This approach is related to recent articles [IBll^, which address 
the numerical solution of singularly perturbed equations of react ion- diffusion type.) In 
a more general numerical-analysis context, we note that sharp estimates for continuous 
Green's functions (or their generalised versions) frequently play a crucial role in a priori 
and a posteriori error analyses [3 [TTl [18] . 

The purpose of the present paper is to establish sharp bounds for the derivatives of 
the Green's function in the Li norm (as they will be used to estimate the error in the 
computed solution in the dual Loo norm [6]). Our estimates will be uniform in the small 
perturbation parameter e in the sense that any dependence on e will be shown explicitly. 



Note also that our estimates will be sharp (in the sense of Theorem 2.4) up to an e- 



independent constant multiplier. We employ the analysis technique used in [H] , which we 
now extend to a three-dimensional problem. Roughly speaking, we freeze the coefficients 
and estimate the corresponding explicit frozen-coefficient Green's function, and then we 
investigate the difference between the original and the frozen-coefficient Green's functions. 
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This procedure is often called the parametrix method. To make this paper more readable, 
we deliberately follow some of the notation and presentation of [8j. 

The paper is organised as follows. In Section [2| the Green's function associated with 



problem ( 1.1 ) is defined and upper bounds for its derivatives are stated in Theorem 2.2 , the 
main result of the paper. The corresponding lower bounds are then given in Theorem 2.4 
In Section [sj we obtain the fundamental solution for a constant-coefficient version of (1.1) 
in the domain f2 = M^. This fundamental solution is bounded in Section HI It is then 
used in Section [5] to construct certain approximations of the frozen-coefficient Green's 
functions for the domains Q = (0, 1) x and Q = (0, 1)^. The difference between 
these approximations and the original variable-coefficient Green's function is estimated 



in Section |6| which completes the proof of Theorem |2.2 

Notation. Throughout the paper, C, as well as c, denotes a generic positive constant that 
may take different values in different formulas, but is independent of the small diffusion 
coefficient e. A subscripted C (e.g., Ci) denotes a positive constant that takes a fixed 
value, and is also independent of e. The usual Sobolev spaces W^'^{D) and Lp{D) on 
any measurable domain D G M.^ are used. The Lp(D)-norm is denoted by IHIp;/) while 
the W"^'P {D)-noTm is denoted by ||-||m,p;Z)- By £c = {xi,X2,X3) we denote an element in 
M^. For an open ball centred at of radius p, we use the notation B{x',p) = {a; G : 
J2k=i23(^k — x'^)^ < p^}. The notation d^^f, ^x„,f A^, is employed for the first- 
and second-order partial derivatives of a function / in variable Xm, and the Laplacian in 
variable x, respectively, while d^^^^f will denote a mixed derivative of /. 

2 Definition of Green's function. Main result 



The Green's function G = G{x;^) associated with (1.1), satisfies, for each fixed x 



LlG{x;$,):=-6A^G + a{$,)d^,G + b{^)G = 6{x-i) foi^eQ, (2.1a) 

G{x-$) = foT^edQ. (2.1b) 

Here is the adjoint differential operator to Lx, and 6{-) is the three-dimensional Dirac 



5-distribution. The unique solution m of (1.1) allows the representation 



u{x) = jjj^G{x-i)f{i)di. (2.2) 

It should be noted that the Green's function G also satisfies, for each fixed ^ G 

LxG{x;i) = -e/\xG-d^,{a{x)G) + b{x)G = 5{x-i) ioixeVl, (2.3a) 

G{x-i) = Q ioixedVt. (2.3b) 

Consequently, the unique solution v of the adjoint problem 

Ll.v{x) = —e AxV + h{x) dx^v + c{x)v = f{x) for £c G ^2, (2.4a) 

v{x) = for a; G dVt, (2.4b) 
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is given by 



v{0 = jjj^G{x-i)f{x)dx. 
We start with a preliminary result for G. 



(2.5) 



Lemma 2.1. Under assumptions (1.2), the Green's function G associated with (1.1) 
satisfies 



If \G{x-i)\di2di,<C, \\G{x- 

^(0,1)2 



1:0 



(2.6) 



'(0,1) 

where G is some positive e-independent constant. 



Proof. The first estimate of (2.6) is given in the proof of [H Theorem 2.10] (see also 
[191 Theorem III. 1.22] and [2J for similar two-dimensional results). The second desired 
estimate follows. □ 

We now state the main result of this paper. 



Theorem 2.2. The Green's function G associated with (1.1), (1.2) in the unit-cube do- 
main Q = (0, 1)'^ satisfies, for all x E Q, the following bounds 

\\d^Mx;-)\\,,n<Gil + \\ne\), (2.7a) 
||%G'(a^; Olli;^^ + \\d.,M^-, ■)hn < Ce-^'\ k = 2,?,. (2.7b) 

Furthermore, for any ball B{x', p) of radius p centered at any x' E Cl, we have 

\\G{x;-)\\,^,,Bi.',p)<Ge-'p, (2.7c) 

while for the ball B{x,p) of radius p centered at x we have 

\\dlG{x;-)\\,,n\B(.^p) < G6-'ln{2 + 6/p), (2.7d) 
WdlGix; ■)\\i;n\B(.,,) < Ge-\\ \ne\ + ln(2 + e/p)), k = 2, 3. (2.7e) 

Here G is some positive e-independent constant. 

We devote the rest of the paper to the proof of this theorem, which will be completed 
in Section |6l 



In view of the solution representation (2.2), Theorem 2.2 yields a number of a priori 



solution estimates for our original problem. E.g., the bounds (2.7a), (2.7b) immediately 
imply the following result. 

Corollary 2.3. Let f{x) = d^,Fi{x) + d^,F2{x) + d,,,F3{x) with Fi, F2, F3 G L^{n). 



Then there exists a unique solution u G Loo(^) of problem (1.1), (1.2), for which we have 
the bound 



\u\ 



<C[(l+|ln£|)||Fi|U;n+£^'/'(||F2||oo;n+||i^3||oo;n)]. 



(2.8) 
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It can be anticipated from an inspection of the bounds for an explicit fundamental 
solution in a constant-coefficient case (see Section |4]) that the upper estimates of Theo- 
rem 



2.2 are sharp. Indeed, one can prove the following result. 



Theorem 2.4 ([9J). Let e G (0,Co] for some sufficiently small positive Cq. Set a{x) := a 



and b{x) := in (1.1). Then the Green's function G associated with this problem in the 
unit cube Vt = (0, Vy satisfies, for all x e [|, |]^, the following lower bounds: 

||%G(a3;-)||i;n > c|ln£|, (2.9a) 
\\d^,Gix; > ce-'^^ k = 2,3. (2.9b) 



Furthermore, for any ball B{x; p) of radius p < |, we have 
wnf^- \\\ ^\cp/e, forp<2e, 

\\G{x;-)\\i,i;nnB{a,;p)> < . ,.1/2 . (2.9c) 

lc(p/£) ' , otherwise., 

\\dlfi{x- ■)||i;f7\i?{x;p) > ce-^ ln(2 + e/p), for p < c,e, (2.9d) 

WdlGix; OlliMBCx;^) > ce'\\n{2 + s/p) + | \ne\) forp<\,k = 2, 3. (2.9e) 

Here c and ci are e-independent positive constants. 



3 Fundamental solution in the constant-coefficient case 



In our analysis, we invoke the observation that constant-coefficient versions of the two 



problems (2.1) and (2.3) that we have for G, can be easily solved explicitly when posed 
in M?. So in this section we shall explicitly solve simplifications of (2.1) and (2.3). To 



get these simplifications, we employ the parametrix method and so freeze the coefficients 



in these problems by replacing a(^) by a[x) in (2.1), and replacing a[x) by a(^) in (2.3) 



and also setting 6 := 0; the frozen-coefficient versions of the operators L*^ and Lx will be 

denoted by L*^ and Lx, respectively. Furthermore, we extend the resulting equations to 
and denote their solutions by g and 'g. So we get 

L\ 9{x; ^) = -e A^g{x; ^) + a{x) %^(x; $,) = 5{x - ^) for ^ G M^ (3.1) 
Lxg{x;^) = -eAxg{x;$,)-a{^)d,,g{x-$) = 6{x-^) for a; G (3.2) 



As X appears in (3.1) as a parameter, so the coefficient a{x) in this equation is considered 
constant and we can solve the problem explicitly. Setting q = ^a{x) for fixed x G (0, 1)^ 
and g{x:,^) = V{x;$,) e'^^^^'^ (see, e.g., [13]), one gets 



-e'A^V + q'V = e e'^^^/^ 6{x - ^) = e e'^^^/" 6{x - $,). 



qxi/e 



As the fundamental solution for the operator —e A^ + q is 
so 



1 6-9''/ = 



Chapter VII], 



Vix; = £e-"i«/" 



1 • 
4^' 



-rq/e 



where r = a/ (xi - ^i)^ + (x2 - 6)^ + (a^s - ^s)^- 
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Finally, for the solution of (3.1) we get 



I Q<l{ii-^i-r)/e 



where q = q{x) = ^a{x). 



A similar argument yields the solution of (3.2) 

_ X Qqi^i-xi-r)/e 



where q = g(^) = |ci(^). 



Let = (^1 - xi)/e, 6 = (^2 - X2)/e, ^3 = (^3 - Xg)/^ and f[^,] = y ^^^^j + ^| + ^1- 

As we shall need bounds for both g and g, it is convenient to represent them via a more 
general function 



9 = 9{x]$]q) :-- 



1 Qq{^i,i.i]-n^i]) 



as 



gix;^) = g{x-i-q) 



1 r ^' 



9{x;i) = g{x;i\q) 



<i=2°-{i) 



(3.3) 



(3.4) 



We use the subindex [xi] in ^ij^i] and r^^^j to highlight their dependence on xi as in many 
places Xi will take different values; but when there is no ambiguity, we shall sometimes 
simply write ^1 and r. 



4 Bounds for the fundamental solution g{x] ^; q) 

Throughout this section we assume that Vl = (0, 1) x M^, but all results remain valid for 
Q = (0, 1)^. Here we derive a number of useful bounds for the fundamental solution g of 
( 3.3[ ) and its derivatives that will be used in Section sj As in ^ and ]) we set q = ^a{x) 
and q = |ci(^), respectively, so we shall also use, for k = 2,3, the differential operators 



(4.1) 



Lemma 4.1. Let x G [—1,1] x and < < q < C . Then for the function 
g = g{x; ^; q) of (3.3) we have the following bounds 



\\9{x; 




l;f7 < C, 




(4.2a) 


\\9ii9{x; 




|l;f7<C(l 


+ llnel). 


(4.2b) 






^ C, 


k = 2,3. 


(4.2c) 


||(CT[^,]%^)(a;; 








(4.2d) 






1 ;n < C*, 


k = 2,3. 


(4.2e) 
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and for any ball B{x'] p) of radius p centered at any x' G [0, 1] x M^, we have 

\\9{x; ■; g) II 1,1 ■nnB(a^';p) < Ce~^p, 
while for the hall B{x] p) of radius p centered at x, we have 
\\dl^9{x;-;q)\\i-n\B{a^;p) < C£"Mn(2 + e/p), 



\\dl9ix; •; q)\\i;n\Bi.;p) < Ce-\\ni2 + e/p) + \ \ne\), 



2, 3. 



Furthermore, one has the bound 

\\dx^9{x; 



1:Q 



<C(l + |lne|), 



and with the differential operators (4.1), one has, for k = 2, 3, 

\\Di:^g{x--;q)\\i.,n + \\D,,g{x;-,q)\U.,n < Ce~^l\ 
||(ef[^,]D5,9^,^)(a;;-;g)||i;Q + IK^n^i] ■; g)||i;n < Ce'^l^ . 



(4.2f) 

(4.2g) 
(4.2h) 

(4.3a) 

(4.3b) 
(4.3c) 



Proof. First, note that Vx9 = —'^^9, so ( 4.3a[) foll ows from (4.2b), (4.3b) follows from 
^A\, ( |4.2c|), while ( |4.3c[ ) follows from (|4.1[), ( 4.2e[ ). Thus it suffices to establish the 
bounds (14.21. 

Throughout this proof, whenever k appears in any relation, it will be understood to 
be valid for A; = 2, 3 (as all the bounds in (4.2) that involve k, are given for both k = 2, 3). 
A calculation shows that the first-order derivatives of g = g{x] ^; q) are given by 



9^1 9 
dq9 



1 



1 



q(r-^i) - 



^1 



r J 



1 



(4.4a) 
(4.4b) 
(4.4c) 



Here we used d^.r = e ^^j/r for j = 1, 2, 3. In a similar manner, but also using 
—£~^^i^j/r^ with i 7^ j, one gets second-order derivatives 



%(e./r) = 

9la9 



1 6 



1 ^ . 



3^1 -r , Ai 



+ 3 



47r£:3 
1 

Atte^ 



r 



- r) + g 
-q&-rf^ — 



Finally, combining c}|^(7 = —d'^^g — d^^g + y t^giS' with (4.4a) and (4.5c) yields 



(4.5a) 
(4.5b) 
(4.5c) 



9l9 



Ane^ 



(4.5d) 
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Now we proceed to estimating the above derivatives of g. Note that d$, = e^d$,, where 
^ eh := e~^{—xi, 1 — xi) X C (—00, 2/e) x M^. Consider the two sub-domains 

Qi := { fi < 1 + |f }, Q2 := { max{ 1, |r } < 6 < 2/s }. 

As r2 C U 1^2 for any xi G [—1, 1], it is convenient to consider integrals over these two 
sub-domains separately. 

(i) Consider ^ e Qi. Then ^1 < 1 -|- |r, so one gets 

s'[{l + r){e~'\g\ + |%^| + \d^,g\ + \d,g\ + \dl,g\) + er\dl^^g\] 

<C^^e-'i^/\ (4.6) 

where we combined e«^i < e«(^+^/2) ^j^h (1 + f + f2 ^ ps) < C'e«?-/4 r|.j^-g immediately 
yields 



Hi 




[(1 + r){e-'\g\ + |%^| + \d^,g\ + \d,g\ + \dl^g\) + sr\dl^,g\\ (e'd^) 

1 

POO 

<C e-''^'Ur<C. (4.7) 
Jo 

Similarly, 



so 




Qi\B(0;p) 



/■oo 

O^li^l + l^l^l] {e^dl)<Ce-' {r-' + r)e-^^'^dr<Ce-'H2 + r')- 
;p) •'p 

(4.8) 

Furthermore, for an arbitrary ball B-^ of radius p in the coordinates ^, we get 

//[ [1^1 + + (^'^^) < C / e-^^/^df < Cmin{p, 1}. (4.9) 

J J Ja,inBp Jo 

(ii) Next consider £ G Jl2- In this sub-domain, it is convenient to rewrite the integrals 
in terms of (^1, t2, h), where 

tk-.^il'^'L so i;'^'dik^dtk and r-i, = ^;^<tl + tl^:t\ (4.10) 

Note that < f < 2 in 1^2 so f - fi = (§ §)/(r fi) > cq^^ where cq := |. 
Consequently e~^^^~^i^ < e~^'^°*^ or 

e-9(5^-6) < CfQ^ where g := ff^ e-^'^"*', (4.11) 
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and 



{i + t + e + e + t^)Q 



{l + t + t^ + t^ + t^) e-"'"'' dt2 dt3 



< C. 



(4.12) 



Using (4. 4), (4. 5) and (4.10) it is straightforward to prove the following bounds for g and 
its derivatives in VL2 

e^\g\<CeQ, 
e'\dlg\<Ce-^i^'[l + t']Q, 



(4.13a) 
(4.13b) 
(4.13c) 



and also 



e\er\d^,g\ + \d,g\)<Ce[l + t^]Q, 

e\eT\d^,^,g\)<cti"^t[l+e]Q, 
e\er\dl^g\)<Ce{e + t')Q, 

e'\dlg\<Ce-^i^\l + e + t')Q. 



(4.13d) 

(4.13e) 

(4.13f) 
(4.13g) 

(4.13h) 



Combining the obtained estimates (4.13) with (4.12) yields 




^2 



[\g\ + e"'\d^,g\ + er\d^M + \d,9\ + e'/hr\dl^^g\ + er\dl^g\ + \dlg\] [e'dl) 



<C j''\e + s"%''^]dl<C. 



Similarly, combining (4.13c) and (4.13e) with (4.12) yields 




^2 



[\diA+e\dlg\] {e'dl)<C j'\'dt<C{l + \\ne\) 



(4.14) 



(4.15) 



Furthermore, by (4.13b), and (4.13e) for an arbitrary ball -Bp of radius p in the coordinates 

we get 




{\g\ + \d^A + {e'dl) < C \e + ff' + fr'^'] < Cp. (4.16) 



To complete the proof, we now recall that Q C Q1UQ2 and combine estimates (4.7) and 
( 4.8[ ) (that involve integration over fii) with (4.14) and (4.15), which yields the desired 
bounds (4.2a)-(4.2e) and (4.2g), (4.2h). To get the latter two bound we also used the 
observation that the ball B{x; p) of radius p in the coordinates $, becomes the ball 5(0; p) 
of radius p = e~^p in the coordinates ^. The remaining assertion (4.2f ) is obtained by 
combining (4.9) with (4.16) and noting that an arbitrary ball B{x'; p) of radius p in the 
coordinates ^ becomes a ball Bp of radius p = e~^p in the coordinates ^. □ 
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Our next result shows that for Xi > 1, one gets stronger bounds for g and its deriva- 
tives. These bounds involve the weight function 



A 



,2g(xi-l)/6 



(4.17) 



and show that, although A is exponentially large in e, this is compensated by the smallness 
of g and its derivatives. 

Lemma 4.2. Let x G [1, 3] xM^ and < |a < g < C. Then for the function g = g{x; ^; q) 
of (3.3) and the weight A of (4.17), one has the following bounds 

\\i[l + er[^,]]\g)ix;-,q)\\,,a < Ce, (4.18a) 
\\{Xd^,g){x; ■;q)\U,n + U\dgg)ix; ■;q)\U,n < C, (4.18b) 

\\{[l + e'/%,]]Xd^,g){x;-;q)\\,,n + e'/^{er[,,]Xdl^^g){^^^ A; = 2,3, 

(4.18c) 

\\r[x,]dg{Xg){x]-]q)\\i.n + \\eri^,]dg{Xd^,g){x] -^q^i-n < C, (4.18d) 
and for any ball B{x'; p) of radius p centered at any x' G [0, 1] x M^, one has 

||(A^)(a;; ■;g)||i,i;nnB(a,';p) < Ce'^p, (4.18e) 
while for the hall B{x; p) of radius p centered at x and k = 1, 2, 3, one has 

(4.18f) 



(4.19a) 
(4.19b) 



II (A dig) {x;-,q)\\, ,n\Bi.,,) < Ce-' ln(2 + e/p) . 



Furthermore, with the differential operators (4.1) and = 2, 3, we have 



\d^^{Xg){x;-;q)\\i,n + ||D^.,(A^)(a;; •; g)||i;n + ||£'^,(A^)(a:;; •; g) ||i ;n < C, 

Ikn^i] D^,^{Xd^,g){x;-;q)\\i.n+ \\er[^^] D^^d-,^{Xg){x; q)\\i.n < Ce"^/^ 



Proof. Throughout this proof, whenever k appears in any relation, it will be understood 
to be valid for k = 2,3 (as all the bounds in (4.18), (4.19) that involve k, are given for 
both A; = 2, 3). 

We shall use the notation A = A{xi) := (xi — l)/e > 0. Then (4.17) becomes 
A = e^'^^. We partially imitate the proof of Lemma 4.1 Again d$, = d$,, but now 
^eh = £-'i-Xi, 1 - xi) X M2 c (-3/e, -A) X So ^1 < -A < immediately yields 



(4.20) 



Consider the sub-domains 



n[:={ \l\<l+^-r, fi<-A}, 

:= { 161 > max{l, |f}, -3/e < < -A }. 

As f2 C n\LiQ,2 for any Xi G [1, 3], we estimate integrals over these two domains separately. 

10 



(i) Let $, e Vt\. Then |^i| < I + \r so, by ( |420| , one has Ae^^^ < e^^^+'^Z^). The first 
inequahty in (4.6) remains vahd, but now we combine it with 



(4.21) 



(which is obtained similarly to the final line in (4.6)). This leads to a version of ( 4.7[ ) that 
involves the weight A: 



A [[l + r){e-M + \d^A + \d^,g\+ e~'\d,g\ + \dl^g\) + er\dl^^g\] [e'dl) < C. 

(4.22) 




In a similar manner, we obtain versions of estimates (4.8) and (4.9), that also involve the 
weight A: 




n[\B{0;p) 



X\dlg\ {6'di)<Ce-Hn{2 + r'), 




n[nBp 



^[\g\ + \d^.g\ + {M) < Cmin{p,l}, 



(4.23) 



(4.24) 



where Bp is an arbitrary ball of radius p in the coordinates Furthermore, (4.22) 



combined with \dg{Xg)\ < \{2A\g\ + \d,g\) and \d,{Xd^,g)\ < X{2A\d^,g\ + \dl^g\) and 
then with A < 2/e yields 




d,{Xg)\+e\d,{Xd^,g)\] {e'di) < C. 



(4.25) 



(ii) Now consider ^ G fig. In this sub-domain (similarly to ^2 in the proof of 
Lemma 4. 1 ) one has iCil < f < 2|^i| and Cot^ < f - < t^, where 4 := 



for k = 2,3, and := ^2 + ^s? (compare with ( |4.10[ )). We also introduce a new barrier Q 

g := A^^e^^^^-I^il) {|eir^e-«^o*'} e'^^^-^^) < Cf Q, (4.26) 



(compare with (4.11); to get the bound for e '^^'^ we used (4.20[)). 



With the new definition (4.26) of Q, the bounds (4.13a)-(4.13c) remain valid in fig 
only with replaced by |^i|. Note that the bounds (4.13d)-(4.13g) are not valid in fig. 



(as they were obtained using r — C^i < t^, which is not the case for < 0). Instead, using 



> 1 and r < 2|^i|, we prove, directly from (4. 4), (4. 5), the following bounds in fig: 



e^\d^,g\<CQ, 
e'\d,g\<CeUQ, 
e\eT\d^,^^g\)<C\W^tQ, 
e\\d,{Xg)\ + e\d,{Xd^,g)\) <CsX -A)+t' + l] Q. 



(4.27a) 
(4.27b) 
(4.27c) 
(4.27d) 
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In particular, to establish (4.27d), we combined dq{Xg) = X[2Ag + dgg] and dq{Xd^^g) = 
X[2Ad^-^g + c^l-^gf?] with the observations that 

(f + iCil) -2A = 2(161 -A) + (f - 161) < 2(161 -A)+t' 

and r~^A <C. 

Next, note that (4.12) is valid with Q replaced by the multiplier { [61''^ 6"*^^°* } from 
the current definition (4.26) of Q. Combining this observation with the bounds (4.13a)- 



(4.13c) and (4.27a)-(4.27c), and also with f < 2|6|, yields 




X [{e''+r)\g\ + l^^^l + {I + s''^r)\d^,g\ + \dqg\ + e''\er\dl^^g\) + e\dlg\\ [e'dl] 



max{A,l} ^ 

<C I [l+eU + + {eUfl' + 161-^] e^'^^^-l^^l) ^6 < C. 

-3/e 



(4.28) 



Similarly, from (4.27d) combined with r < 2|6| < ^, one gets 




n',, 



r[\dq{Xg)\+e\dq{Xd^,g)\] {e'd^) 



max{A,l} ^ 

< C I [(161 -A) + l] e'«(^-l«il) rf6 < C. (4.29) 

-3/e 



Furthermore, by (4.13b), and (4.27a), for an arbitrary ball -Bp of radius p in the coordinates 
?, we get 




max{A,l} 

A [|^?| + \d^,g\ + \d^,g\] {e'di) <C j [l + U'^"] e^^^^-I^^D^ifi 



fignSp J- max{A,l}-p 



< Cp. 



(4.30) 



To complete the proof of (4.18), we now recall that Vt C Vt'^UVt^ and combine estimates 



( |4.22D , ( |4.23D , ( |4.25[ ) (that involve integration over with ( |4.28D , ( |4.29), whi ch yields 
the desired bounds (4.18a)-(4.18d) and the bounds for d'^^g and ^Igf? in (4.18f). To get 



the latter two bounds we also used the observation that the ball B(x;p) of radius p in 
the coordinates ^ becomes the ball 5(0; p) of radius p = e~^p in the coordinates The 
bound for d'^_^ g in (4.18f) follows as d'^_^g = —d^^g — d'^^g d ^^g for ^ x. The remaining 
assertion (4.18e) is obtained by combining (4.24) with (4.30) and noting that an arbitrary 



ball B{x'; p) of radius p in the coordinates ^ becomes a ball Bp of radius p = e ^p in the 



coordinates Thus we have established all the bounds (4.18) 



We now proceed to the proof of the bounds (4.19). Note that Vxg = — V^(7. Combin- 
ing these with (4.18b) and the bounds for ||Ac?^2(y'||i;f2 and ||A 9g3(7||i;Q in (4.18c), yields 



+ \\XD.^^g\\i.n + \\XD^^g\\i.n < C. 
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Now, combining dx^X = 2q6 



and dqX 



2 AX < Ae-^X with dilSal), yields 



\gdcc^X\\i.n + \\gDx^X\ 



1 ;Q 



\gD^,X\U,n<C. 



Consequently, we get (4.19a) 



To estimate erDx^. (A d^^g), note that it involves erdx^, (A d^^ 
we have a bound in (4.18c), and also erdq{Xd^^g), for which we have a bound in (|4.18d|). 



—e9X d'^-^^^g for which 



The desired bounds for erDx^.{Xd^-^g) in (4.19b) follow 



For erD^^dx^{Xg) in (4.19b), a calculation yields e'rD^^,dxi{Xg) = er D^^{Xdx^g) + 
2r^D^^{qXg). The first term is estimated similarly to er^Dxf,{Xd^-^g) in (4.19b). The 
remaining term 7^D^j^{qXg) involves r^d^^{qXg) = qr^Xd^^g, for which we have a bound in 
(4.18c), and also rdg{qXg) = qrdq{Xg) + rXg, for which we have bounds in (4.18d) and 
(4.18a). Consequently (4.19b) is proved. □ 



Lemma 4.3. Under the conditions of Lemma for some positive constant Ci one has 
l|A^(a^; •)ll2,i;[oi]xR^ + \\Dx,{Xg){x- Olli.i^piixM^ < Ce-'=^"/^ A; = 2, 3. (4.31) 



Proof. We imitate the proof of Lemma 



:= n', n [i, < 



4.2 



only now < | or < (| — Xi)/e < 
Thus instead of the sub-domains ^2'^ and fig we now consider Vt'[ and fig defined by 
Xi — \)/e}. Thus in fi'/ (4.21) remains vahd with q > ^a, but now 



r > l/e. Therefore, when we integrate over fi" (instead of fi'^^), the integrals of type 



(4.22), (4.23) become bounded by Ce ^^"^^ for any fixed Ci < ^. Next, when considering 



integrals over fig (instead of fig), note that A — < — 1/£ so the quantity e'^'^^^ l^^'-' in 

Consequently, the integrals of type 

□ 



the definition (4.26) of Q is now bounded by e" 



(4.28) over fig also become bounded by Ce '^i"/'^. 



Remark 4.4. The estimates of Lemmas 4-1 o,nd 4-2 remain valid if we set q := ^a{x) or 



q := |a(^) in g, X, and their derivatives ( after the differentiation is performed) 



5 Approximations G and G for the Green's function G 

We shall use two related cut-off functions ujq and ui defined by 

a;o(t) G C^(0,1), a;o(t) = 1 for t < |, a;o(t) = for t > |; uji{t) := ujq{1 - t), {h.l) 

so Ujn{m) = 1, ujm{l - ui) = and uj'^{t)\^^^^^ = uj'^it)\^^^^^ = for m = 0,1. 

Our purpose in this section is to introduce and estimate frozen-coefficient approxima- 
tions G and G of G. We consider the domain fi = (0, 1) X ]R2 in 

the first part of this 

section, and the domain fi = (0,1)^ in the second part. Note that although G and G will 
be constructed as solution approximations for the frozen-coefficient equations, we shall 
see in Section [6] that they, in fact, provide approximations to the Green's function G for 
our original variable-coefficient problem. 
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5.1 Approximations G and G in the domain = (0, 1) x 

To construct approximations G and G, we employ the method of images with an inclusion 
of the cut-off functions of (5.1). So, using the fundamental solution g of (3.3), we define 

(5.2) 



47re2 



. r[2-xi] r[2+xi] 



a;i(6)j, (5.3a) 
a;o(xi) ^ . (5.3b) 



Note that _^ ^ = and _q ^ = (the former observation follows from r[a.j] = ^j.^.^] 
at = 0, and r^^i] = ^2-2:1] ^^"^ H-^'i] ~ ^2+2:1] ^1 ~ shall see shortly (see 

Lemma 5.1 ) that L^G ~ L^G and LjcG !=a Lj-G; in this sense G and G give approximations 



for G. 

Rewrite the definitions of Q and ^ using the notation 



gid\ ■■= g{d,X2,X3]$]q) 

X± ■- g2g(l±xi)/e^ 

and the observation that 

1 Qii^i.i^ir'^id]) 



I Q<}{il.ld]-r[d]) 



r 



P 



[d] 

-2qx\/e 



Id] 



for d = ±xi, 2 ± Xi. 



They yield 



1^ 



G{x;^;q) = [g[xi] -P9[-xi]] - [A g[2-xi] -pA+^p+^.j] a;i(^ 
§{x;^;q) = [g[x,] - A^^p-^i]] - [p9[-x^] - p A+^p+^-i]] wo(a;i). 



(5.4a) 
(5.4b) 

(5.5) 



(5.6a) 
(5.6b) 



Note that A^ is obtained by replacing xi by 2 ± xi in the definition (4.17) of A. 
In the next lemma, we estimate the functions 

^{x; := LIG - L^G, 4>{x- := L^G - L^G. 



(5.7) 



Lemma 5.1. Let x E Q = (0,1) x M^. Then for the functions and of (5.7), one has 

Uix; ■)||i,i;f^ + WdxJix; OIkn + \\dx,Hx; -ni-^n + Mx; OHuiO < Ge'^^"/^ < G. (5.8) 
One also has 



{x;$,)ki^dn = 0. 



(5.9) 
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Proof, (i) First we prove the desired assertions for 0. By (5.2), throughout this part of 
the proof we set q = ^a{x) 



Recall that g solves the differential equation (3.1) 



with the operator LJ. Comparing the explicit formula for g in (3.4) with the notation 
( |5.4aD implies that L*^[rf] = - d) 5(6 - Xa) 5(6 - X3). So, by (|2 l] ), L^gi^,] = L^G, 
and also L'^g[d] = for (i = — Xi, 2 ± xi and all ^ G Now, by (5.6a), we conclude that 



^^[^1(6)^2] where ^2 := A fi-p-xi] -pX^g[x^+2], and L^^2 = for ^ G 
From these observations, = 2eu[{C,i)d^-^Q2 + [^1^1(6) ~ 2ga;^(6)]^2- The definition 

vanishes at 6 



(5.1) of Ui implies that 



assertion (5.9). Furthermore, we now get 



and for 6 ^ |- 



This implies the desired 



x: 



\i,i;n + \\dx2<Pix; + ||5x30(a;; ■)\\i;n 



<C(\\g2(x; 



12,1 ;[Oi]xl 



+ \\Dx.g2(x; 



+ \\D.^.,g2(x; 



ll,l ;[0|]xR 



Combining this with the bounds (4.31) for the terms )^^g[2±xi] of Q2, and the observation 
that 1-0x2^1 + \Dx3P\ < C\dgp\ < C and d^^^p = 0, /c = 1, 2, 3, yields our assertions for (p in 
(0. 

(ii) Now we prove the desired estimate (5.8) for (p. By ( 5.2[ ), throughout this part of 
the proof we set q = |a(^) > ^a. Comparing the notation (5.4a) with the explicit formula 



for g in ( pR| ), we rewrite ( pOj ) as L^g[x^] = 6{x- So L^g[x^] = L^G, by ( p^. Next, 
for each value d = —xi,2 ± xi respectively set s = — 6,=F(2 — 6)- Now by (3.3), one 

I gg{s-il)/£-5nd] 

hasf[rf] = a/(s - + (6 - X2)2 + (6 - 2:3)7^ so5((a:;;s,6,6;?) = — o ^ • 

Note that Lxg{x] s, 6, 6; = ^{^1 ~ ^{^2 — 6) ^{^3 ~ ^3) and none of our three values 

-g^Cfij^il-Hd])- 



of s is in [0, 1] (i.e. 6{s 



0). Consequently, 



for all X e Vt. 



Comparing ( |5.3b ) and (5.6b), we now conclude that (j) = —-^x [1^0(6)^2] where Q2 '■= 
P9[-x^] -pX^g[2+xi] and LxQ2j= ioi x e n. ^ 

From these observations, = 2eco'Q{xi) dx^Q2 + [^1^0 (^1) + 2q'Wo(^i)]^2- As the defini- 



tion (5.1) of Wo implies that vanishes for Xi < |, we have 



U{x;-)\\i,i;n<G max \\d':M^--)h,i;n. 

a!e[|,l] xIR2 

fc = 0,l 



Here G2 is smooth and has no singularities for xi G [|, 1] (because ?^[2+xi] > H-^i] — 1^ ^ 
for X G [|, 1]). Note that \\dl^g[-x{\\\i,i;n < Ge'"^, and \\d^^{X^g[2+x^])\\i,i;Q < Ge'"^ (these 
two estimates are similar to the ones in Lemmas 4A_ and |4.2| but easier to deduce as 
they are not sharp). We combine these two bounds with \d^-^d^^-^^d^d'^^p\ < G£~'^p = 

Ce^^e"^^^^/*^ for fc, / + m + < 1. As for xi > | we enjoy the bound e'^^^i/"^ < e~5°/^ < 

Ce^e~^"/^, the desired estimate for follows. □ 

Lemma 5.2. Let the function R = i?(x;^) be such that \R\ < C minjeqa,], 1}. The 
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functions G and G of ( 5.2[ ), ( |5.6 ) satisfy 



\\G{x;-)\U,n + \\Gix; 

\\d^,G{x; 

l/2||^ D a2 



|l;n<C^, 

|i;n<C^(l + |ln£|), 
|l;Q<C£-l/^ A; = 2, 3, 
|l;Q<C, fc = 2, 3, 



||(i?%G)(=r;-)||i;f. + ^^/1l(i?9|^^,G)(=r; 

and for any ball B{x'\ p) of radius p centered at any x' G [0, 1] x M^, one has 

\Gix; ■)\i,i-B{x';p)nQ. < Ge^^p, 
while for the ball B{x; p) of radius p centered at x, we have 

\\dlG{x--)\\r,n\Bi.;p) < Gs-'\ni2 + e/p), 

\\dlG{x; ■)\\i;n\Bi.;p) < Ge-\H2 + e/p) + | In^l), A; = 2, 3. 



Furthermore, we have for k = 2, 3 



\i:n 



< Ge-^l\ 



\\d^,G{x--)\\^,^<Ge-^l\ 

f {\\iRdl^^G){x; ■)\\^,n + \\d.M^; Olko) dx, < Ge~^l\ 
Jo 



(5.10a) 
(5.10b) 
(5.10c) 
(5.10d) 



(5.10e) 



(5.10f) 
(5.10g) 



(5.10h) 
(5.10i) 

(5.10j) 



Proof. Throughout the proof, whenever k appears in any relation, it will be understood 
to be valid for A; = 2, 3. 

First, note that r[_xi] > H^i] and r[2±xi] > H^i] for all ^ G therefore 



\R\ < G min{ef[^j], eri_^^], eri2-xi], £^2+^^], l}. 



(5.11) 



9[±xi 



Note also that in view of Rer nark 4.4 all bounds of Lemma 4J^ apply to the components 

4.2 apply to the components \'^g[2±xi] of Q and Q in (5.6). 



and all bounds of Lemma 



Asterisk notation. In some parts of this proof, when discussing derivatives of Q, we 
shall use the notation Q* prefixed by some differential operator, e.g., dx^G*- This will 
mean that the differential operator is applied only to the terms of the type g[d±xi], e.g., 
dx^G* is obtained by replacing each of the four terms g[d±xi] in the definition (5.6a) of Q 
by dx^g[d±xi] respectively. 



1. The first desired estimate (5.10a) follows from the bound (4.2a) for g[±xi] and the 
bound (4.18a) for X'^gi2±xi] combined with \p\ < 1 an d \uJo ^i\ < 1 (in fact, the bound 
for G can obtained by imitating the proof of Lemma 2.1). 



2. Rewrite (5.6a) as 



g = gi- wi(6)^2, where Qi := g[x^] -pg[-x^], G2 ■= \ fi-p-xi] -pX^g[2+x^]- 
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As q = ^a{x) in G (i.e. p and in G do not involve one gets 

Now the desired estimate ( |5.10b ) follows from the bound (4.2b) for d^j^g[±x-t^], the 
bound (4.18b) for d^-^g[2±xi], and the bound (4.18a) for X^g[2±xi]- Similarly, our 
next assertion (5.10c) follows from the bound (4.2c) for d^f,g[±xi], and the bound 
( |4.18c| ) for X^d^^^g[2±x^]- 
The next estimate (5.10d) is deduced using 

\Rd^,G\ < \Rd^,gi\ + c\d^,g;\ + c\g2\, IRdl.Gi < \Rdl.g*\ + c\d^,g*2\. 



Here, in view of (5.11), the term Rd^^gi is estimated using the bound ( |4.2d ) for 
er^±x-i]di,i9[±x^]i while the terms R d^^^^Q * are estimated using the bound (4.2e) for 
£^n±2^i]^|i5t_fi'[±^:i] ^^'^ the_ bound ( |4.18c[ ) for er]2±x^\dl^^^g[2±x{\- The remaining 
terms d^^^g^, ^2 and d^^^gi appear in d^^G and so have been bounded when 

obtaining (|5.10b[), (|5.10c[). 



3. The next assertion (5.10e) is proved similarly to (5.10b) and (5.10c), only using the 
bound (4.2f) for g[±x^] and the bound (4.18e) for X^g[2±x^\- 



4. As g = \a{x) in G, then d'^^G = <9|^^*, m = 1, 2, 3, and the assertions (5.10f) and 
(5.10g) immediately follow from the bounds (4.2g) and (4.2h) for d'^^^g[±x-i^] combined 
with the bounds (4.18f) for X^d'^^g[2±xi] where m = 1, 2, 3. 



5. As g = ^a{x) in G, so using the operator Dx,^ of (4.1), one gets 

dx,G = Dx^[g[x,] - P gi-x^]Y - u^ii^i) [Dx^{\~ gi2^x^]) - P Dx,{X^gi2+x, 

-ldx,a{x) ■ dqp ■ [gi^x,] - uJi{^i)X^ g[2+x-,]] 



where \dqp\ < C by (5.4b) (and we used the previously defined asterisk notation). 
Now, dx^G is estimated using the bound (4.3b) for Dxf.g[±xi] the bound ( |4.19a ) 
for Dx,^{X'^g[2±xi])- For the term g[-xi] in dxf.G we use the bound (4.2a), and for the 
term X^g[2+xi] the bound (4.18a). Consequently, one gets the desired bound (5.10h) 



for Dx,G*. 

To estimate Rdc ^ G, k 



2, 3, a calculation shows that 

^l^xp = iDxAi)[9[xi] ~P9l-x^]Y- wi(^i) [Dx^i\'d^^gi2-xr]) -pDx,{X^d^^gi2+x, 
-ldx,a{x) ■ dgp ■ - uJi{^i)X+d^,g[2+x^]] - ^^(6) dx,G2 

where G2 := 



I g=a(x) /2 ' 



. The assertion ( |5.10h ) for Rd^^^^G is now deduced as fol- 
lows. In view of (5.11 ), we employ the bound (4.3c) for the terms £r[±xi]Dxf,d^j^g[±x^] 
and the bound (4.19b) for the terms £r[2±xi] Dx^iX"^ d^j^g[2±xi])- For the remaining 
terms (that appear in the second line) we use \R\ < C and \dqp\ < C. Then we 
combine the bound (4.2b) for ^£_-^g[-x^] and the bound (4.18b) for d£^g[2+xi\- The 
term dx^G2 is a part of dxf^G, which was estimated above, so for dx^G2 we have the 
same bound as for dx^G in (5.10h). This observation completes the proof of the 



bound for Rdl^^G in ( [slOhl ) 



G in ( 


5.10h) 


(5.10h 


)■ 
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6. We now proceed to estimating derivatives of G, so q = ^a(^) in this part of the 
proof. Let ■= 9[±xi] - ^'^g[2^xi]- Then ( |5.6b[ ), ( |5.4b[ ) imply that G = -poG~, 



where po := ujo{xi)p = ujo{xi)e '^i^^l'^ . Note that 

Combining this with \{—2xi/e)po\ < Ce""^^^/^ and q>\a yields 

|%Po| < C, / (|9,,po| + \DAPo\) dx, < [ (Ce-^e-i^^i/^) dx^ < C. 
Jo Jo 

(5.12) 



Furthermore, we claim that 

\\G-\\i;n<C, \\d,,G^h.,n<C{l + \\ne\), 



\D^,G^\\i.n<Ce-'/^. (5.13) 



Here the first estimate follows from the bounds (4.2a), (4.18a) for the terms g[-xi] 
and X^gi2+xi]- The estimate for dxj^G"^ in ( |5.13[ ) follows from the bound (4.3a) for 



dx^g[±xi] and the bound (4.19a) for dxj^{\^g[2±xi])- Similarly, the estimate for D^^G^ 
in ( |5.13 ) is obtained using the bound (4.3b) for D^^g[±x^] and the bound (4.19a) for 

Dil{X^g[2±xr])- 

Next, a calculation shows that 



% G = D^^G^ - Po D^^ G -D^^po-G' 



dx^G = dx^G^ - Po dx^G - dx^po ■ G~ 



Combining these with (5.12), (5.13) yields (5.10i) and the bound for dx^G in (5.10j) 



To estabhsh the estimate for Rd^^^^G in (5.10j), note that 

dl^^^G = D^^dx.G^ - Po ■ D^^dx.G' - dx.Po ■ D^^G' - d^^po ■ dx,G' - D^^dx.Po ■ G~ 



In view of (5.11), (5.12) and (5.13), it now suffices to show that \RD^J)x]G^\\\.yi < 
Ce'^l'^ . This latter estimate immediately follows from the bound (4.3c) for the terms 
^r{±x^\Di^^dx^g\±x^\ and the bound ( |4.19b[ ) for the terms er^^^i] -0^-5x1 (A=^5'[2±xi]) • 
This completes the proof of (5.10j ). 



□ 



5.2 Approximations for the Green's function G in the domain 

n = (0,1)3 

We now define approximations, denoted by G^ and G®, for the Green's function G in our 



original domain f2 = (0, 1)^. For this, we use the approximations G and G of (5.2), (5.3) 
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for the domain (0,1) x and again employ the method of images with an inclusion of 
the cut-off functions of (5.1) in a two-step process as follows: 



G{x-i) -c^o(6)G'(a^;6,-6,e3) 

Gn{x-i) - Wo(^3)G'n(a;;^i,^2,-6) 

G{x]i) - UQ{x2)G{Xi,-X2,X^]i) 

Ga{x;i) - UJQ{X'i)Ga{xi,X2,-X^]i) 



-6^1(6)^(3;; 6, 2 -6, 6), 
-a;i(e3)G'o(a;;ei,e2,2-e3), (5.14a) 

- UJi{x2) G{Xi, 2 - X2, X3; ^), 

- Wi(x3) Gn(xi, X2, 2 - X3; i). (5.14b) 



Then G, 



w\ 



Ci=o,i 



and Gi 



w\ 



1x1=0,1 

furthermore, by (|5.1l), we have 



(as this is valid for G and G, respectively), and 



€fc=o,i 



and Gri 



lxfc=0,l 



for A; = 2, 3. 



Remark 5.3. Lemmas 5.1 and 5^ of the previous section remain valid ifVt is understood 



as (0, 1)^, and G and G are replaced by G^ and G^, respectively, in the definition (5.7) 
of and (f) and in the lemma statements. 

This is shown by imitating the proofs of these two lemmas. We leave out the details 
and only note that the application of the method of images in the ^2- o,nd ^3- (x2- and 



x^-) directions is relatively straightforward as an inspection of (3.3) shows that in these 
directions, the fundamental solution g is symmetric and exponentially decaying away from 
the singular point. 

As G^ and G^ in the domain Vt = (0, 1)^ enjoy the same properties as G and G in the 
domain (0, 1) X R^, we shall sometimes skip the subscript ® when there is no ambiguity. 



6 Proof of Theorem for Q = (0, 
(general variable-coefficient case) 



We are now ready to establish our main result. Theorem |2.2| for the original variable- 



coefficient problem (1.1) in the domain fl = (0, 1)^. In Section[5| we have already obtained 
various bounds for the approximations G® and G® of G in = (0, 1)^. So now we consider 
the two functions 

v{x;0 := [G-Ggp](x;C), v{x;^) = [G - G^]{x; 

Throughout this section, we shall skip the subscript ® as we always deal with the domain 
fi = (0,1)3. 

Note that, by ( |5.7[ ), we have LxV = La,[G — G] = [L^ — Lx]G — (p, and similarly 
L^v = L|[G — G] = [L^ — L^]G — (f). Consequently, the functions v and v are solutions of 
the following problems: 



Lxv{x; ^) = h(x] ^) for x G Vl, v{x; ^) = for a? G dVl, 
Llv{x; ^) = h{x; ^) for ^ e n, v{x; ^) = for ^ G dn. 



(6.1a) 
(6.1b) 
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Here the right-hand sides are given by 



where 



h{x; ^) := d,, {R G}{x; ^ - b{x) G{x; ^) - <P{x- ^) 
Hx; ^) := ^) - 6(0 Gix; ^) - ^{x; 



i?(a;; ^) := a(a;) — a(^), so |-R| < C min{£:f[^.j], 1}. 



(6.2a) 
(6.2b) 



(6.3) 



Applying the solution representation formulas (2.2) and ( 2.5[ ) to problems (6.1a) and 
(6.1b), respectively, one gets 



v{x;$,) 
v{x;$,) 





G{x; s) h{s] ^) ds, 
G{s] ^) h(x] s) ds. 



(6.4a) 
(6.4b) 



We now proceed to the proof of Theorem 2.2 



Proof. Throughout the proof, whenever k appears in any relation, it will be understood 
to be valid for k = 2,3. 



(i) First we establish (2.7b). Note that, the bounds (S.lOi) and (5.10h) for d^f.G and 
^x^^G, respectively, it suffi ces to show that ||%u(a3; + \\^x^^v{x■, ■)\\i-n < Ge^'^l'^. 

Applying (9^^, to (6.4a) and to (6.4b), we arrive at 



%^(»^;^) = 
d^^v{x-X) = 

From this, a calculation shows that 





G{x; s) d^^^h{s; ^) ds, 
G{s;$,)da:^h{x;s)ds. 



\\d^^^v{x;-)\\i.Q < i sup // \G{x;s)\ds2ds3] ■ sup ||%^(s; ■)||i;n c^Si , 
^sie{o,i) J Jm.^ ' Jo (S2,S3) 

\\dxk^ix;-)\\i.n < ( sup ||G'(s; Olli;^-) ■ \\d^Mx'r)\\i;n- 
^ sen ^ 



So, in view of (2.6), to prove (2.7b), it remains to show that 

\\d^Ji{x--)\\^,^<Ge-^l\ 



/ sup \di^Ji{x\-)\i.p,dxx<Ce 

^0 {x-i,x-i)&?- 



These two bounds follow from the definitions (6.2a), (6.3) of h and /i, which imply that 



\d^,h{x-i)\ < \Rdl^^G\+G{\d,,G\ + |%G|) + 1%^ 
\d^M^;^)\ < \Rdl^G\+G{\d^,G\ + \d^,G\) + \d^,^ 
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combined with the bounds (5.8) for 0, 0, the bounds (S.lOi), (S.lOj) for G and the bounds 
(5.10b), (S.lOh) for G. Thus we have shown (2.7b). 



(ii) Next we proceed to obtaining the assertions (2.7a), (2.7d) and (2.7e). We claim 
that to get these bounds, it suffices to show that 

^2 



V := max sup \\df v(x; 
k=2,3 ^en" 



sup 



\i-n<C{l + \\ne\+eV), 



snp\\dlvix;-)\U,n<CE-'il + EV). 



(6.5a) 
(6.5b) 
(6.5c) 



Indeed, there is a sufficiently small constant c=k such that for e < c*, combining the bounds 
(6.5a), (6.5b), one gets W < C(l + |ln£:|), which is identical with (2.7a). Then (6.5a) 
implies that V < Ce~^, which, combined with (5.10g), yields (2.7e). Finally, V < Ge^^ 
combined with (6.5c) and then (5.10f) yields (2.7d). 

In the simpler non-singularly-perturbed case of £ > c*, by imitating part (i) of this 
proof, one obtains W < Gi, where Gi depends on c*. Combining this bound with (6.5a) 



and (6.5c), we again get (2.7a), (2.7d) and (2.7e 



We shall obtain (6.5a) in part (iii) and (6.5b) with (6.5c) in part (iv) below 
(iii) To get (6.5a), it suffices to set k 



2 and consider V := d'^^v (as d'^^v is estimated 



similarly). The problem (6.1b) for v implies that 

LlV{x; ^) = H{x; ^) for ^ e n, V{x; ^) = for ^ G dn. 



The homogeneous boundary conditions di v 



follow from v\ 



(6.6) 

in ( |6.6 ) for m = 1, 3 immediately 



=0,1 



0. The homogeneous boundary conditions on the boundary edges 



,^2 = 0, 1 are obtained as follows. As v\ 
again m = 1, 3. Combining this with h 



6=0,1 

1 6=0,1 



so d^j^vl 



6=0,1 



™ 16=0,1 



0, where 
5.3i we 



6=0,1 = 0- 
d^^hi + h2 



(for which, in view of Remark 

used (5.9)) and the differential equation for u at ^2 = 0, 1, one finally gets ff^^v 

For the right-hand side H in (6.6), a calculation shows that H = H{x] ^) = 
with hi = hi{x]$,) and h2 = h2{x;$,) defined by 

h ■■= d^Ji - 2%a(^) ■ := a|a(^) • - 2%6(^) ■ %t; - 9|6(^) ■ 

with k = 2. Here we used d'l^[ad^^v] = ad^y + 2d^^ad'l_^^J) + d'^^ad^-^v = ad^^y + 
[2 d^^ad^^v] — 9|^a d^^v and [hv] = bV + 2 d^J)d^^v + d^Jjv. (Note that H is under- 
stood in the sense of distributions; see Remark 6.1 below.) 



Now, applying the solution representation formula (2.5) to problem (6.6), and then 
integrating the term with hi by parts, yields 



V{x-^) 




[-ds^G^s; $,) hi{x; s) + G{s; $,) h2{x; s)] ds, 



(for the validity of the above integration by parts we again refer to Remark |6.1 ). As 
(2.7b) implies sup^g^ \\ds2G{s; ■)\\i;n < Ge^^l"^ , while (2.6) implies sup^g^ 11^(5; ■)|| < C, 



imitating the argument used in part (i) of this proof yields 



|l;n 



\\V{x--)h,^<G{e-^l^\hi{x- 



|l;n 



\h2{x\ 
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So to get our assertion (6.5a), it remains to show that 



(6.7) 



To check this latter bound, note that \hi\ + |/i2| < C{\d^i,h\ + \d^-^^v\ + \d^f.v\ + \v\) with 
k = 2. Note also that 



\v{x--)\\,^,,n<C{e 



-1/2 



W) + \\G{x; 



1.1 ■,n, 



where we employed v = G — G and then the bounds (2.6), (2.7b) and the definition (6.5b) 
of W for G. Combining these two observations with 



< \Rdl^^G\ + + + |G|) + 



k 



(where we used (|6.2b), (6.3)), and then with the bounds (5.10a)-(|5.10d) for G, and the 



bound (5.8) for 0, one gets the required estimate (6.7). Thus (6.5a) is established. 

(iv) To prove (6.5b) and (6.5c), rewrite the problem (6.1b) as a two-point boundary- 
value problem in ^i, in which x, ^2 and ^3 appear as parameters, as follows 



[sdl + a(^)%] v{x; = l{x- for ^ G (0, 1), v{x- = 0, 



where 



l{x- i) := h{x- i) + e [dlv{x- i) + dlv{x- i)] - b{^) v{x; ^). 



(6. 



(6.9) 



Consequently, one can represent v via the Green's function F = r(,^i, ,^2, ^3; ■§) of the 
one-dimensional operator [—£d'^_^ + a(^)(9gj. Note that F, for any fixed ^3 and s, 
satisfies the equation [—£d^_^ + a(^)5gjr(^; s) = 6{C,i — s) and the boundary conditions 



16=0,1 



0. Note also that 



[]d^J{^;s)\d^,<2a-' 
Jo 



(6.10) 



[H Lemma 2.3]; see also [13 (1.1.18)], [Hi (3.10b) and Section 3.4.1.1]. 
The solution representation for v via F is given by 



v{x]^) 



T{^;s)h{x;s,^2,^3)ds. 



Applying d^-^ to this representation yields 



sup 



\d^,msMi]- h{x 



1;^ 



In view of (6.10), we now have ||i;f7 < 2a ^\\h\\i.^Q,. Note that the differential equation 
(6.8) for V implies that £1191^^111;^ < Cdl^^j-uHi;!^ + ||/i||i;Q). So, furthermore, we get 

\\d^M\i;n + e\\dlv\\i-n < C||^||i;f^. 
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As G = v + G and we have the bound (5.10b) for d^^G^ to obtain the desired bounds (6.5b) 



and (6.5c), it remains to show that ||/i(cc; < C + eV. Furthermore, the definitions 



(6.9) of h and (6.5a) of V, imply that it now suffices to prove the two estimates 

IKX^;-)lli;n<c, \\Kx--)\\^,^<G. 



(6.11) 



(2.7d) and (2.7e) 



The first of them follows from v = G — G combined with (2.6) and ( |5.10a ). The second is 
obtained from the definition (6.2b) of h using (5.10h) for ||/29^jG||i.f7, (5.10a) for 
and (5.8) for ||0||i;f7. This completes the proof of (6.5b) and ( |6.5c ), and thus of (2.7a), 



(v) We now focus on the remaining assertion (2.7c), again rewrite the problem (6.1b) 



as 



\dn 



0, 



where 



hoix; ^) := h{x; ^) - a(^) d^M^; ^) + [1 - v{x; ^). 



(6.12) 



We shall represent v via the Green's function \l/ of the two-dimensional self-adjoint 
operator [— eA^ + 1]. Note that \E' = \Ef(s;^), for any fixed s, satisfies the equation 
[— eA^ + l]^'(s;^) = 6{$, — s), and also the boundary conditions ^(*j^)|^ggj^ = 0. Fur- 
thermore, for any ball B{x'; p) of radius p centred at any x' , we cite the estimate [3l 
(3.5b)] 



|*(s; ■)\i,i;B{x';p)nn < Ge V- 
The solution representation for v via \E' is given by 



(6.13) 



v{x]i) 




\l'(s; ^) hQ^X] s) ds. 



Applying 9t , m = 1, 2, 3 to this representation yields 



x\-)\i.i;B(^';p)nn < (sup |^(s; ■)|i,i;B(^/;p)nn) • \\ho{x] 



(6.14) 



To estimate ||/io||i;n; recall that it was shown in part (iv) of this proof that H^^iwHi ;q < 
2a;~"'^||/;.||i;n and ||/i(a;; ■)||i;f7 < G + eV, and in part (ii) that V < Ge~^. Consequently 
ll'^^i'^^lli;^ < C. Combining this with (6.12) and (6.11) yields ||/io||i;f7 ^ C. In view 



of (6.14) and (6.13), we now get \v\i^i-B{x';p)nQ. ^ Ce ^p, which, combined with (5.10e) 
immediately gives the final desired bound (2.7c). 



Remark 6.1. Note that the term 9? h in H, where k 



□ 



2, 3, has such a singularity at 



^ = £c that it is not absolutely integrahle in VL. So H and the differential equation in (6.6) 
are understood in the sense of distributions fTR Chapters 1, 3]. In particular d'^^h is a 
generalised ^f.- derivative of the regular function d^Ji. 
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